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We present a three-flavour fit to the recent v^ — > v e and v^ — V v^ T2K oscillation 

data with different models for the neutrino-nucleus cross section. We show that, even 

for a limited statistics, the allowed regions and best fit points in the (#13, <$cp) and 

(#23; A- m atm) planes are affected if, instead of using the Fermi Gas model to describe 

£Nj | the quasielastic cross section, we employ a model including the multinucleon emission 

k* ■ channel. 
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1 Introduction 

Recently the T2K collaboration has released data in both z/ M — > v e appearance [1] and 
v il -^ v ii, disappearance [2] modes; in the first case, six events passed all the selection criteria, 
implying (under the assumption of a normal ordering of the neutrino mass eigenstates): 

sin 2 (2^ 13 )T2X = 0.11, (1) 

with the CP phase Sep undetermined. In the disappearance channel, the 31 events collected 
by T2K are fitted with: 

(sin 2 2# 23 )t2^ = 0.98 \Am 2 atm \ T2K = 2.65 ■ 10~ 3 eV 2 . 

The aim of this work is to reanalyse the T2K data to assess the impact of different models for 
the ^-nucleus cross sections on the determination of oscillation parameters. This work can 
be considered as a generalization of Ref. [3], where the impact of different modelizations of 
quasielatic cross sections in the low-gamma beta-beam regime was analyzed. In the present 
case we consider two different models involving not only quasielastic but also pion production 
and inclusive cross sections. On one hand, we choose a model as similar as possible to the 
one used by the T2K collaboration. They simulate the neutrino-nucleus interaction using 
the NEUT Monte Carlo Generator [4J. Even if we do not know the details of the last tunings 
performed by the collaboration to take into account for the recent measurements of K2K 
[5JE], MiniBooNE [3, [8] and SciBooNE 0,[1Q], we treat the several exclusive channels using 
the same models implemented in NEUT. As a consequence, we consider the Fermi Gas 
[TT] for the quasielastic channel and the Rein and Sehgal model [T2j for pion production. 
The second model considered in our analysis is the one of Martini, Ericson, Chanfray and 
Marteau [13], in the following called "MECM model". It is based on the nuclear response 
functions calculated in random phase approximation and allows an unified treatment of 
the quasielastic, the multinucleon emission channel and the coherent and incoherent pion 
production. The agreement with the experimental data in the pion production channels [6j 
7, 9J has been proved. Nevertheless the main feature of this MECM model is the treatment of 
the multinucleon emission channel in connection with the quasielastic. In fact, as suggested in 
[T3l fT4"] . the inclusion of this channel in the quasielastic cross section is a possible explanation 
of the MiniBooNE quasielastic total cross section |8], apparently too large with respect to 
many theoretical predictions [15] employing the standard value of the axial mass. Since 
the MiniBooNE experiment, as well as many others involving Cherenkov detectors, defines a 
"quasielastic" event as the one in which only a final charged lepton is detected, the ejection of 
a single nucleon (a genuine quasielastic event) is only one possibility, and one must in addition 
consider events involving a correlated nucleon pair from which the partner nucleon is also 
ejected. This leads to the excitation of 2 particle-2 hole (2p-2h) states; 3p-3h excitations 
are also possible. Nowadays other models [TH \TT\ HH [19] have included the multinucleon 
contribution in the computation of the cross sections relevant for the MiniBooNE quasielastic 
kinematics, improving the agreement with the experimental data. For a brief review see for 
example [2D]- Recently, it has been shown [21] that the MECM model can also reproduce the 
MiniBooNE flux averaged double differential cross section [8] which is a directly measured 
quantity and hence free from the model-dependent uncertainties in the neutrino energy 
reconstruction, and the total inclusive cross section [20J (also employed by T2K as described 
below) measured by SciBooNE [10]. In the following we will use the cross sections obtained 
in the two different approaches described above in several exclusive channels (quasielastic 
and pion production), as well as in the inclusive one, for both charged current (CC) and 



neutral current (NC) interactions on carbon and oxygen (the targets used in near and far 
T2K detectors, respectively) and for two neutrino flavours z/ M and v e . Although all exclusive 
channels are involved in the analysis, we will refer to the first model as "the Fermi Gas 
model" and to the second approach as "the MECM model" . 

In order to perform our comparison among the above-mentioned models, we first need to 
correctly normalize the Fermi gas to the T2K event rates, at both near (ND) and far (FD) 
detectors; we use the following algorithm: 

1- normalization of the cross section with the v^ inclusive CC at the ND; according to pQ, 
we have to reproduce 1529 v^ inclusive events, collected using 2.9 x 10 19 POT, in the 
energy range [0 — 5] GeV, with an active detector mass of 1529 Kg Q| at a distance of 
280 m from the v source and half a year of data taking (Run 1). Notice that only the 
muon neutrino cross sections can be correctly normalized; we assume that the same 
normalization also applies for the v e cross section, although they could differ at the /i 
production threshold (in any case away from the peak of the neutrino flux); 

2- computation of the expected events (and energy distributions) at the far detector in 
the appropriate two-parameter plane ((sin 2 2# 13 , 5c p) for appearance and (#23, Am 2 4m ) 
for disappearance); 

3- normalization to the T2K spectral distributions. 

Step #3 is needed to get rid of the experimental efficiencies applied by the T2K collaboration 
to the signal and background events. This means that the bin contents of our simulated 
distributions (obtained at point #2) are corrected by coefficients, generally of 0(1) that we 
consider as a detector property, and then not further modified. For a different model, we 
repeat step #1 and then go to step #2, using the same normalization coefficients extracted 
in step #3 with the Fermi gas. We make use of the GloBES [22] and MonteCUBES [23] 
softwares for the computation of event rates (and related \ 2 functions) expected at the T2K 
ND and FD detectors. The fluxes of v^ v e and their CP-conjugate counterparts predicted at 
the FD in absence of oscillations have been extracted directly from Fig.l of [lj, whereas the 
v^ flux at the ND has been obtained from [2]. Such fluxes (the relevant ones summarized in 
FigJT]) are given for 10 21 POT. As already stressed, for the relevant cross sections we assumed 




Figure 1: Fluxes at near (v^ only) and far (v^ and v e ) detectors. 



: We thank Scott Oser for providing such a number to us. 



that the T2K collaboration uses some "sophisticated" version of the Fermi gas model |llj . 
In Figj2]we show the inclusive and QE cross sections in the FG model (dashed lines) and in 
the MECM model (solid line) used in our simulation, after having correctly normalized the 
inclusive cross sections to the event rate at the ND. Especially for the MECM model, this 
procedure involves a degree of extrapolation of the inclusive cross sections towards neutrino 
energies beyond the validity of model itself. However, neutrino fluxes above 0(1) GeV drop 
very fast and we checked that different kind of extrapolations do not alter our conclusions. 
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Figure 2: Inclusive (thin lines) and QE (thick lines) v^ CC cross sections on oxygen in the 
FG model (dashed lines) and in the MECM model (solid line) after the normalization of the 
inclusive cross sections to the event rate at the ND. 

The important feature here is that, even after the normalization procedure, the MECM 
CCQE cross section is still larger than the FG predictions, in the energy range relevant for 
appearance studies. This is due to the inclusion of the multinucleon component and will be 
the main reason of the differences between the results obtained in the two models. Note on 
the contrary that the inclusive cross sections are not really different. 

2 The appearance channel 

The U/j, —¥ v e transition probability is particularly suitable for extracting information on 
013 and Sep] at the T2K energies (E u ) and baseline (L), one can expand the full 3-flavour 
probability up to second order in the small parameters #13, A 12 /A 13 and A 12 L, with A^- = 
Am 2 ,/4£, 



where 



s 2 3 sin 2 2#i3 sin 2 (A atm L) + c 2 3 sin 2 26 12 sin 2 (A so ; L) 
J COS (S C P + A atm L) {A sol L) sin (2 A atm L) , 



J = C13 sin 2#i2 sin 2# 2 3 sin 2#i 3 , S23 = sin 6- 



23 



(2) 



(3) 

We clearly see that CP violating effects are encoded in the interference term proportional to 
the product of the solar mass splitting and the baseline, implying a scarce dependence of this 
facility on Sep when only the v^ — > v e channel (and the current luminosity) is considered. 



2.1 Extracting the T2K data 



Events in the far detector (obtained with 2.9 x 10 POT) are u e CCQE from z/ M — y v e 
oscillation, with main backgrounds given by v e contamination in the beam and neutral 
current events with a misidentified 7r°. The experimental data have been grouped in 5 
reconstructed-energy bins, from to 1.25 GeV and they are summarized in TabJT] The 
expectations for signal and backgrounds have been computed by the T2K collaboration from 
Monte Carlo simulations, for fixed value of the oscillation parameters, namely sin 2 2# 12 = 
0.8794, sin 2 2^3 = 0.1,sin 2 2# 23 = 1 and Am 2 sol = 7.5 x 10~ 5 eV 2 , Am 2 atm = +2.4 x lO^eV" 2 . 
In order to normalize our event rates to the T2K Monte Carlo expectations, we extracted 
these numbers from Fig. 5 of [1] and reported them in TabJU 
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0.204 
0.148 


0.783 
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0.08 


0.192 

0.113 

0.04 


4.171 
0.662 
0.565 



Table 1: Expected event rates for sin 2 2# 13 =0.1. 



Notice that we used the central bin energy as a reference value for the neutrino energy 
in a given bin; this could be different from the reconstructed neutrino energies used by 
the T2K collaboration. To mimic possible uncertainties associated to the neutrino energy 
reconstruction, we apply an energy smearing function to distribute the rates in the various 
energy bins. Other choiches, more related to microscopical calculations [25| |26| 127] are 
also possible. In particular, an analysis of the validity of the approximation contained in 
the identification of the reconstructed neutrino energy via a two-body kinematics with the 
real neutrino energy has been done in [27], where the MECM model has been employed. 
The role of several nuclear effects such as Pauli blocking, Fermi motion, RPA correlations 
and multinucleon component has been studied in details. This analysis was performed, 
among others, considering the T2K conditions at near and far detectors, paying a particular 
attention to the v^ — » v e appearance mode. The ratios among our computation and the T2K 
data {energy dependent efficiencies) are summarized in Tabf2J This procedure (corresponding 
to step #3 of the previous paragraph) allows us to take into account all the detection 
efficiencies to different neutrino flavours in the Super Kamiokande detector. Once computed, 
these corrective factors are used in the simulations done with a different cross section, since 
we assume here that they are features of the detector and not of the neutrino interactions. 
We see that for v e ^ — y v e transitions these numbers are just 0(1) coefficients, which makes 
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Table 2: Efficiencies computed after normalizing the event rates at the values for 
sm 2 28 13 =0.1. 



us confident that the normalization procedure correctly accounts for the main experimental 



features. The same is not true for the NC events which, however, have not been normalized 
to the ND as for the CC interactions. As a check, we also computed the expected events 
for sin 2 26 l i3=0, obtaining 0.1 v^ — > v e events and 0.72 v e — > v e events (and the same neutral 
current rate), in good agreement with the T2K expectations pp. 



2.2 Fit to the data 

Equipped with these results, we performed a \ 2 analysis to reproduce the confidence level 
regions in the (sin 2 2#i3, <5cp)-plane shown in Fig. 6 of [lj. Contrary to what has been done in 
the official T2K paper, we make a complete three-neutrino analysis of the experimental data, 
marginalizing over all parameters not shown in the confidence regions. As external input 
errors, we used 3% on 6* 12 and Am^, 8% on #23 and 6% on Am 2 tm . We use a constant energy 
resolution function a(E u ) = 0.085 and, for simplicity, we adopt a 7% normalization error 
for the signal and 30% for the backgrounds. We also used energy calibration errors fixed to 
10 -4 for the signal and 5 • 10 -2 for the backgrounds; normalization and energy calibration 
errors take into account the impact of systematic errors in the x 2 computation. 
Assuming a normal hierarchy spectrum, the best-fit point from the fit procedure is (obvi- 
ously) : 



sin 2 (2#i 3 ) = 0.108 5 C p = 0.04 



(4) 



with x 
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vain 



1.69; the related contour plot is shown in Figj3j Compared to the official 
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Figure 3: The 



sin 2 (20 n ) 

and 90% C.L. regions for (sin 2 29 13 , <5 CP ) in the FG model. 



release, the plot is in quite good agreement, although the allowed values of 9 13 around 
maximal CP violation 5c p = tt/2 are a bit larger (this is the effect of including the errors of 
the atmospheric parameters into the fit procedure). 

We now apply the same procedure to determine 6^3 using the MECM cross sections 
described in [13]. In doing that, we normalize the cross sections to the ND events and then 
compute the number of oscillated events (and related backgrounds), to be compared with 
the experimental T2K data. We assume that the efficiencies reported in Tabj2] are exactly 
the same, since they are a property of the SK detectors and then independent on the cross 
section model. With these assumptions, we get the following number of expected rates for 
sin 2 2013=0.1: 
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Table 3: Total rates for sin 2 29i3=0.1 in the MECM model. 



It is clear that larger rates need smaller #13 to reproduce the data (the effect of the CP 
phase 5 is negligible with such a statistics). The best fit point is: 
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sin 2 (2#i 3 ) = 0.073 5 C p = , (5) 

1.53, and the contour plot is shown in FigJU We can appreciate a substantial 
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Figure 4: The 68% and 90% C.L. regions for (sin 2 29 13 , #cp) for the MECM model. Star 
indicates the best fit point. 

improvement in the determination of the reactor angle, whose largest value is 0.24, to be 
compared with 0.35 obtained with the Fermi Gas. To make a more direct comparison on 



#13 between the FG and MECM results, in Figj5]we show the x 2 



X,r, 



function, computed 



marginalizing over all other oscillation parameters (including Sep)- At la, we get: 
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sm 2 29[ 3 G 
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They are clearly compatible although, as expected, 6^ ECM < 9 



FG 

13 • 



3 The disappearance channel 

We extend the previous analysis to include the first disappearance v^ — > v a data [2]. In the 
two-flavour limit, (the one where both # 13 and Am^ are vanishing) the v n — > v n probability 
reads [28]: 



PU, 



1 - sin 2 26*23 sin 2 (A atm L) 



(7) 
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Figure 5: y 2 as a function of 9\3 for the MECM model (solid line) and FG (dashed line). 

Effects related to #13 are clearly sub- dominant, so that this channel is particularly useful 
to extract information on the atmospheric parameters. The T2K collaboration collected 31 
data events, grouped in 13 energy bins, as one can see from Fig.3 of [2J. The sample extend 
up to 6 GeV and it is mainly given by z/^CCQE, z^CC non-QE, v e CC and NC. Differently 
from the appearance channel, we cannot normalize their energy distribution to the channel- 
by-channel T2K Monte Carlo expectation since, as far as we know, such information has not 
been released. We can only normalize our FG cross section to the total rates shown in Tab. I 
of [2J, which amounts to 17.3, 9.2, 1.8 and <0.1 events for z/ M CCQE, z^CC non-QE, NC 
and v e CC, respectively. Such numbers refer to sin 2 (26*23) = 1.0 and |Am 2 fm | = 2.4 x 1CT 3 
eV 2 , with all other neutrino mixing parameters vanishing. For the sake of completeness, we 
summarize in TabJH the T2K data as well as the energy distributions of the z/^CCQE and 
z^CC non-QE as obtained using the MECM cross sections. In evaluating such numbers, 
we assume a variable bin size, centered in the neutrino energy corresponding to the T2K 
data. In our fit procedure we have assumed a total normalization for the NC as given in |2J, 
but with appropriate energy distributions according to the FG and MECM cross sections. 
We have also adopted a conservative 15% normalization error and energy calibration error 
at the level of 10~ 3 for both signal and background. The results of our fit procedure can 
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Table 4: T2K events and bin distributions for the v^CCQE and v^ CC non-QE rates in the 
MECM model. 



be appreciated in Figj6j where we show the 90% CL for the Fermi Gas (dashed line) and 
the MECM model (solid line), in the case of normal hierarchy. We plot the 2 degrees 
of freedom (dof) confidence levels in the (#23, Am 2 im ) (left panel) and (sin 2 26*23, Am 2 im ) 
(right panel, to facilitate the comparison with the official T2K results) planes. Again, the 
plots have been obtained marginalizing over the not shown parameters (a full three-flavour 
analysis); we considered a 50% error on sin 2 2^13 (with best fit at sin 2 26*13 = 0.0059) and 
5cp undetermined. We obtained: 
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Figure 6: 90% contour levels for the MECM model (solid line) and FG (dashed line), in the 
(#23, Am 2 tm ) (left panel) and (sin 2 2# 2 3, Am 2 tm ) (right panel) planes. Star indicates the best 
fit obtained in the MECM model. 



FG : sin 2 2# 9 , > 0.86 



MECM : sin 2 2# 23 > 0.91 
with best fit points: 

FG : sin 2 2# 23 = 0.99 (47.9° 



2.22 • KT 3 < Am 2 tm (eV 2 ) < 2.90 • 10~ 3 
2.31 • 10~ 3 < Ami ( e V 2 ) < 2.93 • 10~ 3 



(8) 



Am ltm 



2.56- 10" 3 eV 2 



(9) 



MECM : sin 2 2# 23 = 1.00 (45.0°) Am 2 tm = 2.62 • 10~ 3 eV 2 



Some comments are in order; first of all, we observe that, for both models, the best fit point 
is different from the T2K ones, which is 

(sin 2 29 23 ) T2K = 0.98 \Am 2 atm \ T2K = 2.65 ■ 10~ 3 eV 2 ; 

this is somehow obvious since we normalized our events to the MC predictions obtained 
for a different set of atmospheric parameters. The MECM cross section gives a better 
determination of both #23 and Am 2 tm , mainly due to the larger statistics than the FG; at 
the same time, the disappearance probability in Eq. (JTJ), for negligible solar mass difference 
and reactor angle, is smaller if the atmospheric mass difference is larger, for fixed sin 2 2# 23 . 
This is what happens here, where a smaller P(v^ — > v^) (and then a larger Am 2 4m ) is needed 
in the MECM model to partially compensate for the larger cross section. 



4 Future perspectives 



The statistics used by the T2K collaboration to make the disappearance study is only a 
2% of the rates expected at the end of the experiment. It makes sense to ask how the 
previous results would modify if the accumulated statistics would be larger than the current 
one. We limit ourselves to consider a number of events with the same energy distribution 
as the experimental ones but bin contents larger by factor of 10, in both appearance and 
disappearance channels. In the analysis of the appearance channel, the (weak) information on 
#13 contained in the disappearance sample should not be neglected (as we did previously); at 
the same time, the dependence on the atmospheric parameters from the appearance sample 
cannot in principle be neglected when studying the disappearance data. For this reason, we 
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prefer to combine both v^ — > v e and ^ — \ v^ oscillation data, and study the sensitivity to 
the reactor and atmospheric parameters as we did in the previous sections, marginalizing 
over the parameters not expressly shown. Notice that such an approach would not give any 
additional information on the mixing parameters if adopted with the current T2K statistics: 
in fact, we see from Figj6]that the uncertainties on #23 and Am 2 im obtained from the T2K 
data are larger than the adopted external errors on these parameters in the appearance 
channel, so that adding the v^ — > v^ data will not improve the sensitivity to #13; similarly, 
the dependence on the reactor angle in Piy^ — > v^) is sub-leading and the impact of the 
disappearance channel in the appearance measurement is also negligible. We stress that 
extracting information on the mixing parameters combining appearance and disappearance 
channels is also mandatory to avoid some inconsistencies emerged in the official T2K fits, 
where |Ara 2 tm | is fixed to 2.4 x 10~ 3 eV 2 in the appearance analysis whereas the best fit 
point obtained from the disappearance data is 2.6 x 10~ 3 eV 2 . The results of our procedure 
are shown in FigJTJ where we display the 90% CL in the (sin 2 2#i3, <5cp)-plane (left panel) 
and (#23, Am 2 im )-plane (right panel) for the MECM (solid line) and FG (dashed one) models 
in the case the current T2K statistics is increased by a factor of 10. The minimum of the \ 2 
still gets reasonable values: we obtain Xmin ~ 20 in the appearance analysis and Xmin ~ ^0 
in disappearance. In both panels we can appreciate the effects of the increased statistics, as 
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Figure 7: 90% CL for 2 dof in the (sin 2 20 13 , 5cp)-plane (left panel) and (6*23, Am 2 4m ) -plane 
(right panel) for the MECM model (solid line) and FG (dashed one) in the case the cur- 
rent T2K statistics is increased by a factor of 10. Stars indicate the best fit values of the 
parameters as obtained in the MECM model. 



expected: for a given model of cross section, the allowed regions are strongly restricted with 
respect to the current situation. The best fit values for Sep are somehow different in the two 
models (Sep ~ and Sep ~ 144° for the FG and MECM models, respectively), although 
statistically non very significant. Such intervals for sin 2 2# 13 , #23 and Am 2 tm are summarized 
in Tab]5] (Sep is obviously still unconstrained). We have checked that, if we only use the 
appearance channel to extract # 13 , the results are sligthly different: although the best fit 
value is practically indistinguishable from the one quoted in Tabj5l the confidence regions 
are a bit larger, with significant overlap with the above mentioned analysis. To see stronger 
effects due to the #13 dependence in the v^ — \ v^ transition, we need a more accurate spectral 
information [29J . Similar conclusions can also be drawn for the disappearance channel: with 
only a factor of 10 more statistics and no appearance contribution, the best fit for the 
atmospheric parameters remain almost the same whereas the 90% CL region for #23 shows 
a smaller lower limit (from 40.1° to 39.8°) in the FG model. 

Finally, we observe that such an increased statistics is necessary to make marginally 



sin^ 20 



13 



^231 



Am2 <m (10- 
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FG [0.041-0.211] (0.105) [40.1-51.3] (47.6) [2.45-2.67] (2.56) 
MECM [0.023-0.154] (0.092) [41.1-49.9] (45.4) [2.49-2.67] (2.60) 

Table 5: 90% intervals for sin 2 2#i3, #23 and Am 2 im; for the MECM and FG models in the 
case the current T2K statistics is increased by a factor of 10. In parenthesis, the best fit 
points. 



incompatible the FG and MECM sin 2 2# 13 results, see FigJHl obtained marginalizing over 
5c p also. In fact, at la we get: 
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Figure 8: x 2 — Xmin as a function of sin 2 2^13 for the MECM model (solid line) and FG 
(dashed line) in the case the event rates are increased by a factor of 10. 
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5 The inverted hierarchy case 

For the sake of completeness, we have repeated the same computations as above under the 
hypothesis that the neutrino mass spectrum is of inverted type (IH). With the current T2K 
statistics, we cannot appreciate huge differences in the results obtained using the two different 
models for the cross section. Then, we limit ourselves here to the case where the stastistics 
is larger by a factor of 10. Our results are summarized in Figj9] and Tab]6j Comparing 
the left panel of Figj9] with the corresponding one in FigJTl we clearly see that an inverted 
spectrum prefers larger values for #13, in both models. The best fit of the CP phases is 
different among the two mass orderings but not really significant. In the atmospheric plane, 
right panel of Figj9j we observe the same pattern as in the normal hierarchy case, that is 
the MECM tends to give a better resolution for both Am 2 im and #23 than the FG model. 
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Figure 9: 90% CL for 2 dof in the (sin 2 26 ) i 3 , 5cp)-plane (left panel) and (6*23, Am^ m ) -plane 
(right panel) for the MECM model (solid line) and FG (dashed one) in the case the current 
T2K statistics is increased by a factor of 10 and the "data" are fitted with the inverted 
hierarchy. Stars indicate the best fit values of the parameters as obtained in the MECM 
model. 
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'29 



13 
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|Am 2 t J(10-W 



FG [0.049-0.241] (0.164) [40.0-51.3] (47.6) [2.34-2.55] (2.44) 
MECM [0.026-0.181] (0.102) [41.1-49.8] (45.4) [2.37-2.56] (2.47) 

Table 6: 90% intervals for sin 2 2#i3, #23 and Am 2 tm , for the MECM and FG models in the 



case the current T2K statistics is increased by a factor of 10 and the "data' 
the inverted hierarchy. In parenthesis, the best fit points. 



are fitted with 



6 Conclusions 

In this paper we have studied the impact of using different models for the neutrino-nucleus 
cross section in the determination of the #13,23 mixing angles and the atmospheric mass dif- 
ference Am 2 im using the recent T2K data, for both appearance and disappearance channels. 
Although the statistics is not large enough to draw definite conclusions, we have seen that a 
more refined treatments of nuclear effects in neutrino interactions can have some impact in 
the achievable precision on the mixing parameters. In particular, the MECM model predicts 
a large CCQE cross section, compared to the FG model, which results in a small #13 needed 
to fit the data in the v^ — > v e channel. At the same time, a larger Am 2 4m is required to 
fit the data in the v^ disappearance channel, since a smaller disappearance probability is 
needed to compensate for the larger cross sections. Interestingly enough, with 10 times more 
statistics the two models tend to give substantial different results in terms of best fit points 
and parameter uncertainties (of course, better than before) but their predictions are still 
compatible to each other. 
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